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ABSTRACT 

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Cluster 
Lensing And Supernova survey with Hubble (CLASH) multi-cycle treasury programs with the Hubble 
Space Telescope (HST) have provided new opportunities to probe the rate of core-collapse supernovae 
(CCSNe) at high redshift, now extending to z « 2.5. Here we use a sample of approximately 44 CCSNe 
to determine volumetric rates, Rcc, in six redshift bins in the range 0.1 < z < 2.5. Together with 
rates from our previous HST program, and rates from the literature, we trace a more complete history 
of Rcc{z), with Rcc = 0.72 ± 0.06 yr“^ Mpc“^ 10“"* at z < 0.08, and increasing to 3.7ti/ yr“^ 

Mpc“^ 10“^ /lyo to z « 2.0. The statistical precision in each bin is several factors better than than 
the systematic error, with significant contributions from host extinction, and average peak absolute 
magnitudes of the assumed luminosity functions for CCSN types. Assuming negligible time delays 
from stellar formation to explosion, we find these composite CCSN rates to be in excellent agreement 
with cosmic star formation rate density (SFRs) derived largely from dust-corrected rest-frame UV 
emission, with a scaling factor of fc = 0.0091 ± 0.0017 Mq^, and inconsistent (to > 95% confidence) 
with SFRs from IR luminous galaxies, or with SFR models that include simple evolution in the initial 
mass function over time. This scaling factor is expected if the fraction of the IMF contributing to 
CCSN progenitors is in the 8 to 50 Mq range. It is not supportive, however, of an upper mass limit 
for progenitors at < 20 Mq . 

Subject headings: supernovae: general, surveys 


1. INTRODUCTION 

Core-collapse supernovae (CCSNe) are the explosive 
end products of massive stars. Despite having similar 
explosion mechanisms, these events have a wide range 
of explosion energies due to their broad mass ranges. 
They are much less useful as cosmological probes, and 
as such, for nearly two decades, have arguably been a bi¬ 
product of large surveys for type la supernovae (SNe la). 
Nonetheless, CCSNe are very much interesting in their 
own right. They are far better physically understood 
than their SN la cousins, and most connections to their 
progenitor stars are far better established. These super¬ 
novae are critical to understanding the formation of dust 
and the chemical evolution in galaxies and the intergalac- 
tic medium. Moreover, as they stem from short-lived 
stars, CCSN rates virtually trace the rate of instanta¬ 
neous star formation, providing an independent tracer of 
star formation in galaxies less sensitive to the ambigu¬ 
ity in extinction corrections that plague far-ultraviolet 
measures, and free from UV-to-lR light reprocessing as¬ 
sumptions (e.g., contributions from “warm” vs. “cool” 
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components) inherent in far-infrared measures (cf. Ken- 


nicutt & Evans 20121. Similarly, tracking CCSN rates 


over cosmological distances gives a compelling and less 
biased vi ew of the cosmic star-form ation rate density his¬ 
tory (cf. Madau & Dickinson||2014). 

Until recently, there have not been sufficient samples at 
significant redshifts to track the CCSN rate history with 
any precision, largely as SN la surveys tend to selectively 
bias against CCSNe as they are, o n average, a few mag¬ 
nitudes less luminous than SNe la ( Li et al.||2011 ). Still, 
there have been a small number of complete rate mea¬ 
sures at z > 0.1, notably from: the So uthern Interme- 


2008) 

(z) = 

et al. 

]2(109| 

et al. 



(z) = 0.3), the Subaru Deep Field (Craur 


(z) = 0.66), the Stoc kholm VIMOS Super- 
iurvey ( Melinder et al. |20l2 (z) = 0.39 and 0.73) 


nova 

and the GOODS HST SJN survey (jDahlen et al. 2012 

(z) = 0.39, 0.73, and 1.1 1). _ 

The CCSN rates from Dahlen et al.| (|2012l D12 here¬ 
after) stand out as the first to extend this rate history to 
z > 1. The D12 sample was collected using the HST Ad¬ 
vanced Camera for Surveys (ACS), and used an entirely 
self-consistent analysis to show unambiguous evolution 
in CCSN rate history, or Rcc{z), extending to z = 1.3. 
The D12 results allowed the first viable comparison to 
consensus cosmic star-formation rate density, or '(/’(z), 
partly derived from rest-trame UV observations of the 
galaxies in the same Great Observatories Origins Deep 
Survey (GOODS) fields. While consistent with what 
would be expected from i/'(z), the D12 Rcc{z) measures 
were dominated by large statistical errors, and nearly 
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equally large systematic biases, the largest of which is 
attributable to the obscuration (and loss) of events due 
to line of sight extinction in the SN host. 

The two extragalactic programs of the recent multi¬ 
cycle treasury program with HST, the Cosmic Assembly 
Near-infrared Deep Extragalactic Legacy Survey (CAN- 
DELS) and the Cluster Lensing And Supernova survey 
with Hubble (CLASH), provided new opportunities to 
probe SN rates in high-z galaxies, both utilizing the Wide 
Field Camera 3 (WFC3) IR channel, in programs de¬ 
signed to achieve well-resolved (^ 0".l resolution) survey 
images, to 5 (t sensitivities of ~ 25.5 mag out to 1.6/xm. 
These have a distinct advantage over the D12 ACS survey 
in that they are less sensitive to extinction in host galax¬ 
ies at z < 1, and more sensitive to events out to z = 2.5. 
Both programs have made some noteworthy discoveries, 
i ncluding GSDlOP ri dRodney et al ||2012|), UDSlOW il 
(Jones et al. ] 2013[ ) and ClND12Col ( [Rodney et al.|2C)14[ )- 
three ot the highest redshift SNe la known; ECSll'i'yl, 
GSD12Hum, and GSD12Qua- three of the highest red- 
shift SNe (of any type) to date ( Rodney et al.|2014 ); and 
CL012Car, CLN12Did, and CLAll'lib- three SlNe test¬ 


ing galaxy cluster lens mass models (Patel et al. 2013). 
The analyses on the ra tes of SNe la i r om th ese s urveys 


have been p resented in Rodney et al. 
et al.|(|2014|). 


(2014) and Graur 


in this paper we present Rcciz) from the SNe discov¬ 
ered in the CANDELS and CLASH programs, extending 
CCSN rates for the first time to z = 2.5, and greatly 
reducing the statistical errors well below the systematic 
errors at all lower redshifts. In Section [2] we review the 
observational details of these surveys. In Section we 
describe the selection of approximately 44 of the 94 dis¬ 
covered events we identify as CCSNe, and the calculation 
of the control-time rates from the sample. In Section 
we present the results from this survey, compare them to 
other rates in the literature and compare the combined 
CCSNe rate history to the star formation rate density 
histories as derived from UV and IR studies. Finally, we 
summarize these results in Section!^ Throughout most 
of this paper (except where indicate in the error analy¬ 
sis) we assume Hq = 70 km s“^ Mpc“^, = 0.27, and 


Da = 0.73, c onsistent with the WMAP7 cosmology (Ko- 
matsu et al.||2011 1. 


2. OBSERVATIONAL DETAILS 
2.1. CANDELS 

The CANDELS program surveyed five well-studied 
deep fields: GOODS-S, GOODS-N, COSMOS, UDS, and 
EGS, over a three-year period, each with different total 
integration times necessary to satisfy “dee p” and “wide” 
comp onents of the near-IR galaxy survey (Grogin et al. 
2011). To accommodate searches for supernovae, the 
observations for each field were broken up into semi- 
regularly cadenced visits (or epochs), with a total of two 
visits to each of the “wide” fields and 10 to each of the 
“deep” fields, separated by ^ 50 days to match SN la 
rise times at z « 1.5. 

The exact details of the SN survey are described 

R14 hereafter), but here we 


Rodney et al. (2014 


provide a review of the survey design. Each single-epoch 
set of exposures consisted of at least four WFG3/IR de¬ 
tector exposures: two in the FlOOVE-passband (77), and 
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Fig. 1.— Effect of search area (in arcmin^), cadence (in days), 
and E160VF-passband sensitivity (magnitudes) on the simulated 
number of events at three redshifts. Area has the greatest effect 
on survey yield, while cadence and depth are secondary factors of 
the same magnitude. 


two in F125IT-passband (J), each to a 5cr limiting Vega 
magnitude of miim{H) = 25.4 and mumiJ) = 25.8, re¬ 
spectively. These passbands were the effective search fil¬ 
ters for the SN survey, requiring that events be present 
at least in both 77-band exposures, with a modest re¬ 
quirement for visibility in the J-band exposures as there 
is the possibility that z > 2.5 SNe la could drop out of 
the bluer passband. 

Volumetric SN rate calculations require simulations of 
each survey field, where in which area, cadence, and 
depth are critical factors. These simulations can be¬ 
come intensive if each pair of repeated pointings in a 
mosaic field are considered individually. Generally sur¬ 
veys have a built-in uniformity convenient to considering 
a few large mosaics over many repeated visits. The CAN- 
DELS fields, however, are somewhat less uniform, being 
composed of nine repeated subareas of very different ca¬ 
dences. 

One could, as we have done in RI4, use an effective 
area for each subarea, each with a specific average ca¬ 
dence. However, better precision can be obtained using 
cadence-area groups that better reflect the full range in 
image combinations. For example, Figure [l] shows the 
effect of varying area, cadence, and depth on the rela¬ 
tive expected yield of a mock survey, where area is the 
principal component, and cadence and depth have sim¬ 
ilar secondary but not negligible effects. By design, the 
CANDELS exposure times within a visit are all the same 
so there is no variance in the survey sensitivity. But 
selecting an average representative cadence that is 50% 
higher or lower than the actual cadence for a field can 
lead to a 50% error in the relative number of expected 
events, if survey area remains unchanged. 

For this analysis we collect pointings of similar ca¬ 
dence, regardless of target field, into 30 effective tar¬ 
get cadence-areas of a virtual, two-epoch survey. These 
cadence-areas are shown graphically in Figure and are 
broken down by field in the Appendix. Treating the 
survey in this way has the added benefit of removing 
a source of confusion in control-time rates (see § for 
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Fig. 2. — Diagram of the 30 effective areas in CANDELS (green 
stems), the 24 in CLASH WFC3/IR (red stems), and the 28 in 
CLASH ACS WFC (blue stems) with similar cadences. The areas 
and their samples for CANDELS, CLASH WFC3/IR, and CLASH 
ACS are each treated separately in the rates analysis. 

rolling surveys, where long-term events can be sustained 
over several epochs of the survey, essentially causing a 
double-counting problem in the tracking of rising and 
declining events. 


2.2. CLASH 


Ove r the same three-yea r period, the CLASH pro¬ 
gram ( [Postman et al.|[2M2 ) targeted 25 galaxy clusters, 
with observations spread over up to 16 epochs to enable 


spread over up to 16 epochs 
supernova discovery. The observational strategy allowed 
for pairs of simultaneous observations with WFC3 and 
with ACS, with one focused on the target cluster, and 
the other on one of four parallel positions, constrained 
by the roll angle of HST. The ~ 6' separation between 
instruments provided some assurance the deep paral¬ 
lel fields were not significantly influenced by the grav¬ 
itational leasing from the target cluster, nor contami¬ 
nated by target cluster members, and therefore essen¬ 
tially equivalent to the “open” fields provided by CAN- 

DELS. __ 

As detailed in (Graur et al.||2014 G14 hereafter), the 
CLASH parallel fields were essentially surveys of the red¬ 
dest filters in each detector-set, each visit of which con¬ 
sisting of 2 exposures in WFC3-iL, and 4 to 6 exposures 
in the ACS F850LP-passband {Z), to a 5a limiting Vega 
magnitude of = 25.4 and mumiZ) = 25.2, for 

up to 4 visits per parallel field. As with CANDELS, 
simulating the total of 132 H and 148 Z-band visits as 
independent template-search pairs would have been too 
computationally expensive, so instead we collected each 
into similar 24 and 28 cadence-area groups for H and 
Z, respectively, as shown in Figure and detailed in the 
Appendix. 


2.3. Image Processing, SN Detection, and 
Classification 

SN candidates were discovered in real time, with fast 
FTP delivery of raw data within a few hours of exe¬ 
cution. Our processing pipeline created searchable dif¬ 
ference images within an hour of delivery, the prod¬ 


ucts of which were examined by an experienced search 
team to select candidates for coordinated followup, usu¬ 
ally within 24 hours of the exposures. The processing 
pipeline has several elements similar to the calibration 
pipeline used in the OTFR processing by OPUS, pro¬ 
viding up-to-date flat field calibrations, up-the-ramp IR 
sampling with cosmic-ray flagging , and utilizing Mul- 
tiDrizzle (Koekemoer et al. 2002|) to provide the final 
sampling to an astrometric grid. Lhe few additional cal¬ 
ibrations include improved intra-exposure source-based 
relative alignment, and modest correction for persistence 
from exposures within our own program, both of which 
reduce false detections in the search images. To construct 
searchable subtraction images, all previous IR imaging in 
regions that overlap a given incoming epoch were drizzled 
to the same astrometric grid, and each subtracted from 
the incoming calibrated data. As there is little variation 
in the HST point spread function with WFC3 and ACS, 
direct subtractions were sufficient without the need to 
apply image convolution techniques. 

SN candidates were identified by human review, and 
tracked over the regular cadences of the surveys. Addi¬ 
tional HST photometric target-of-opportunity observa¬ 
tions were made for many events to gather the most com¬ 
plete multi-wavelength light curves possible for as many 
candidates as possible, with priority given to candidates 
likely to be SNe la at z > f. 

Redshift constraints are crucial for developing efficient 
followup strategies, and the benefit to searching these 
well-studied fields has been the copious additional in¬ 
formation on host galaxies, including existing catalogs 
)ic and photometric redshifts (cf. jPahlenl 


et al. 

2013| 

moer 

et al.| 


and multi-wavelength photometry (Koeke¬ 


moer et al.l 20111 from which to derive new photometric 


All SNe in this sample were photometrically identified 
using STARDUST based on their observed light curves, 
and best estimates of their redshifts, photometric or spec- 
troscopic. As discussed in R14, STARDUST uses the 
SNANA (Kessler et al.||2009) light curves and SED tem¬ 
plates to simulate events of specific types at the redshift 
of the candidate. Then, through a Bayesian analysis 
these simulated SNe are compared against the observed 
data to determine the most likely type and condition 
(e.g., peak luminosity and extinction) for each candi¬ 
date event, along with an appropriate range of certainty. 
While STARDUST is principally designed for identify¬ 
ing SNe la, the Bayesian analysis also provides a mea¬ 
sure of a candidate’s likelihood of being some other SN 
type, normalizing such that the sum of probabilities for 
all SN types is unity. The tool only weighs these possibil¬ 
ities, and assumes the likelihood of other identifications 
(e.g., AGN, variable stars, or noise) is negligible or com¬ 
pletely ruled out by other information, such as a lack of 
event-host separation, lack of periodicity, or extremely 
low S/N. It is therefore useful to use the Pia comple¬ 
ment, or 1 — P/q, as the probability of a given event 
being a CCSN. 


3. DETERMINING VOLUMETRIC RATES 

We determined the Rcc in this analysis using the 
control-time method, which has some similarities to the 
monte carlo simulations used by SNANA and other tools, 
e.g., D12 and R14, but differs in how the effective time 
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of the survey is calculated. The volumetric rate of CC- 
SNe in any redshift interval, represented by a mean or 
effective redshift, z, is expressed as: 


Rcc{z) 


N{z) 

tc{z)AV{z)’ 


( 1 ) 


where N{z) is the number of observed CCSNe in the 
redshift interval, and tc{z) is the effective control time in 
the same redshift interval, or similarly the control time 
at the effective redshift for the interval. These two terms 
are discussed further in the subsequent sections. AV (z) 
is the comoving volume sampled in the redshift interval 
(i.e., zi < z < Z 2 ), in the solid angle, fl, of each survey 
field, such that AV{z) = n/47r[y(z2) — V(zi)]. 


3.1. Observed Number of Events 

While in principle it seems easy to count discoveries 
and sort them by type and redshift, these crucial bits of 
information are often ambiguous. Few events in our sam¬ 
ple have spectroscopic redshifts, determined with cer¬ 
tainty from strong features in their hosts’ spectra, and 
fewer still have spectroscopic SN classifications. We rely 
on photometric redshift techniques for many events, and 
the STARDUST photometric classifier for most events, 
to determine probabilistically the most likely scenario 
of type and redshift. We account for this impreciseness 
in sample in our rate calculation using an approach de¬ 
scribed in R14, determining the number of events in any 
redshift bin, N{z), as the product of the photometric 
redshift probability distribution function, or P(z|D), and 
the “CCSN” classification probability, or 1 — Pja, for 
each event, summed for all events, and integrated over 
the redshift interval. In this way: 



N 

r 1 

N{z) 

=/s 

-21 :! i 

Pi(z'|D) X (1-Pla.i) 


Events with spectroscopic redshifts have P(z|D) that are 
delta-functions, and those with spectroscopic classifica¬ 
tions, SN la or CCSN, have Pia = 1 or 0, respectively. 

To allow a more direct comparison to previous rates, 
we calculate our CCSNe rates in six redshift bins cen¬ 
tered at redshifts 0.3, 0.7, 1.1, 1.5, 1.9 and 2.3, each with 
a bin width of Az = 0.4. Both the cumulative redshift 
distribution, and the binned cumulative number distri¬ 
bution used in the valuation of the rates are shown in 
Figures and 1^ 


3.2. The Control Time 

The control time method addresses the probability of 
detecting a SN of a given type, redshift, peak luminosity, 
age, and extinction expressed as an amount of time in 
which the said event would be visible, given the depth, 
area, and cadence of the survey. It starts in very much 
the same way that monte carlo methods do, utilizing li¬ 
braries of light curves and spectra to simulate SNe with 
assumed conditions, under the constraints of the assumed 
distributions of each of these conditions. The monte carlo 
approach, however, samples from these distributions ide¬ 
ally placing simulated events in real search data to test 
their recovery, often in real time. The control time ap¬ 
proach is to first determine survey sensitivities (perhaps 



Fig. 3. — The effective number of CCSNe from the CANDELS 
survey as a function of redshift. Left ordinate: the solid line shows 
the redshift distribution function for all events, determined from 
P(z|D) and CCSN type probabilities, and shown in log space to 
accentuate the redshift distribution. Right ordinate: the green 
region shows the binned cumulative number distribution, or N[z). 



Fig. 4.— Same as Figure]^ the effective number of CCSNe from 
the CLASH survey as a function of redshift. The blue and orange 
lines are for the ACS and WFC3/IR parallel fields, respectively, 
scaled by values in legend to highlight their relative contributions 
to the total redshift distribution function (in black), and binned 
cumulative number distribution (in green). 


in real time), and then iteratively test the recoverability 
of all model scenarios against this sensitivity threshold, 
yielding more accurate but computationally intensive re¬ 
sults. 

The control time (tc) itself is essentially the amount of 
time our survey is sensitive to CCSNe of a given type at 
a specified redshift. In an analytical approach, we de¬ 
termine the likelihood of detecting a model CCSN type 
considered at a light-curve age of t by differencing the 
model flux at this age from the model flux at an ear¬ 
lier age, given by a time-dilation corrected cadence. We 
compare this flux difference to our survey efflciency in 
detecting objects of similar differenced flux, resulting in 
a likelihood, e, for detecting the CCSN model at that 
light-curve age. We integrate over all plausible light- 
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curve ages, from explosion to several hundred days after 
explosion in the rest frame of the modeled event, giving 
the total effective time that SN model would be visible 
within the survey, or, 

OO 

tciP) = j e{l3;t)dt. (3) 

t=o 

Here, /3 holds the adjustable parameters of the model 
(SN type and redshift) and the survey (cadence, hlter, 
and sensitivity). 

The survey efficiency was determined from empirical 
tests of fake SNe in real data (see R14 and G14 for de¬ 
tails). We express it here as 

e(m) = {1-I-exp[0.4(m — TOc)]}“^, (4) 

where m is the magnitude corresponding to the differ¬ 
ence in flux between the detection and template epochs, 
and TOc is the magnitude cutoff, where transients with 
m > rric were recovered < 50% of the time. This value 
is m^{H) = 25.4 for CANDELS and the WFC3/IR com¬ 
ponent of CLASH, and rndZ) = 25.2 for the ACS/WFC 
component of CLASH. 

Even within a given type, CCSNe have a range of 
peak luminosities (characterized by H-band peak abso¬ 
lute magnitude, Mb), and they are likely to experience 
some extinction internal to their host galaxy, which we 
characterize by the extinction parameter in the H-band, 
or Ay- To account for each of these possibilities (dis¬ 
cussed in further detail in later sections), we evaluate 
both as changes to the observed flux, and determine the 
expectation of tc using probability distribution functions, 
P{Mb) and P{Av), respectively, such that 

OO OO OO 

tciP) = J J J e{j3-,t,MB,Ay) P{Mb) 

0 0 —OO 

X _P(Ay) dM^Q dAy dtj (5) 

thereby dehning a control time in the rest frame for an 
entire CCSN type at that redshift in our survey. As the 
details of both surveys have been discussed in Section!^ 
and more extensively in R14 and G14, we will focus the 
following discussion on the selected CCSN models and 
their dependencies. 


3.3. CCSN Types 

CCSNe have very different light curve shapes, peak lu¬ 
minosities, and relative rates of occurrence resulting in 
different visibilities in our survey. For individual assess¬ 
ment, we separate the CCSNe into five main types: Type 
H Plateau (HP), Type H Linear (HL), Type H with nar¬ 
row emission lines (Hn), and stripped envelope Types, 
Ib and Ic. We adopt mean absolute magnitudes and dis¬ 
persions determined from low- z samples from the Asi - 
ago Supernova Catalog (AHC, Richardson et al. 2014), 
reproduced here in Table and use them as Caussian- 
distributed luminosity functions in P[Mb)- 

The ASC sample relied heavily on discoveries as re¬ 
ported by the community, and therefore is likely subject 
to discovery biases of from varied (and perhaps unknown) 
survey cadences. In co ntrast, the Lic k Observatory Su¬ 
pernova Search (LOSS, Li et al.|[2011 ) had a strategy of 


TABLE 1 

Subtypes of Core Collapse Supernovae 


Type 

Peak Mb^ 

Dispersion 

Fraction'^ 

IIP 

-16.80 ± 0.37 

0.97 

52.4% 

IIL 

-17.98 ± 0.34 

0.90 

7.3% 

Iln 

-18.62 ±0.32 

1.48 

6.4% 

Ib 

-17.54 ± 0.33 

0.94 

6.9% 

Ic 

-16.67 ± 0.40 

1.04 

17.6% 

Faint CCSNe 

> -15 


9.4% 


Richardson et al. 


1 2014[|, where uncertainty is the 


^ From 

standard error m tne meam Dispersion is the standard 
deviatio n of the sample . 

From Li et al. (20111 


systematically repeated visits to the same target galaxies, 
which suggests a greater completeness in polling CCSN 
diversity than ASC. However admittedly, ASC may com¬ 
pensate by nature of having a larger volume-limted sam¬ 
ple, to 100 Mpc rather than 60 Mpc. We nonetheless 
adopt the volume-limited, Malmquist-bias corrected frac¬ 
tions of CCSNe for each type from LOSS as the relative 
frequency from each contributing to the total number of 
expected events. Those fractions are also shown in Ta¬ 
ble [D 

The SNANA library includes some well studied CC¬ 
SNe, each with well-sampled ugriz light curves and 
temporally matched optical (3000 A to 1 pm) spectra. 
We use these light curves and spectra to build mean 
templates for each of the CCSN types, ultimately con¬ 
structed from 28 SNe HP, 2 SNe Hn, 11 SNe Ib, and 
8 SNe Ic from the Sloan Digital Sky Survey and the 
Carnegie Supernova Project, and a synthetic SN HL 
model (cf. the SNANA handbook). We normalize these 
light curves to the average peak Mb for each class from 
Table m 

For simplicity, we assume the luminosity functions 
of each CCSN subtype are Gaussian [N{MB,cr)\, with 
means and dispersions from Table and treat these 
as probability distributions for Equation such that 
P{Mb) = N{Mb,(j)- We also modify the ^served dis¬ 
persion to account for an additional scatter due to grav¬ 


itational leasing, which we discuss in Section 3.6 
As discussed in D12, there are populations of rare CC¬ 
SNe that are not well accounted for in the literature. 
For example, there is a population of SN 1987A-hke and 
other less luminous CCSNe [Mb > —15) that do not 
have sufficient numbers for assessment in LOSS or ASC 
samples. Estimates on the fr action faint CC SNe pro¬ 
duced range to as high as 30% ( Li et al.||2011 1, however, 
there is not yet a complete sample to adequately char¬ 
acterize the luminosity function of these subtypes in the 
low-z universe. But, by being intrinsically faint, their 
contribution to our CCSN rates is expected to be lim¬ 
ited to lower redshifts, perhaps contributing only to our 
lowest redshift bin. Conservatively, we estimate their 
contribution at ~ 10%. 


There ar e also classes of super-lum inous events (Gal- 
Yam 2012 |Heger & Woosley 2002), whose origins re¬ 
main uncTeaFTuF'areTTkeljTTied'^Eo extremely massive 
stars, which have the potential to add to every redshift 
bin. However, as they are exceedingly rare occurrences 
at low-z (perhaps < 2%), we chose not to account for 
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these in our expected samples. We note, however, that 
some super-luminous events may stem from population 
III stars, implying that their number could increase with 
redshift. As this possibility is only weakly explored in 
the literature, we leave this as a caveat to this study. 


3.4. K-Corrections 

The template light curves are generally not well- 
matched to the part of the rest-frame spectral energy 
distributions (SEDs) covered by the H and Z-passbands 
of the surveys. We apply k-corrections to the light curve 
templates to simulate the flux in this synthetic rest-frame 
passband. We calculate the k-c orrections following the 
prescription of Kim et al. (1996), 


k{t) = - 2.5 


-log 


log 


/ g(A[l-fz])f(A,t) dX 


JSy{X)FiX,t)dX 

f H{X)Z{X)dX'' 
fSy(X)Z(X)dX^ 


( 6 ) 


where F{X,t) is the spectral energy distribution of the 
SN at the time of observation, H{X) is the transmission 
of the survey passband, F[{X[l + z]) is a synthetic survey 
filter blue-shifted to the rest frame of the SN, and Sy{X) 
is a rest-frame template passband {ugri or z) that pro¬ 
vides the closest match to the synth etic filter. Z(X) is 
the Vega spectral energy distribution ( Colina et aHl996 1 
that defines the zero-magnitude base tor the magnitude 
system. 

3.5. Extinction Correction 

For most of this analysis, we consider CCSN hosts as 
normal galaxies well-represented by the sample of CCSN 
hosts seen in the local universe. At low-z, the observed 
distributio n of host extinction. Ay, for CCSNe is ex - 
ponential (Hamuy & Pinto 2002 Schmidt et al. 1994), 
which is supported by dust modeling of host ga laxy ex¬ 
tinction effects on the observed r ates of CCSNe (Hatano 
et al.|[l9^ Riello fc Patat|[^005[). This ext inction distri¬ 


bution IS corrected using the 
tion law (with Ry = 4.05) 


Calzetti et al. (2000) extinc- 
Tom the rest-frame V-band 


to the effective wavelength of their synthetic rest-frame 
filter covered by the survey passbands. We evaluate this 
distribution as the extinction likelihood function, where 
P{Ay) oc exp(—Ay Ay) in the control time calculation of 
Equation!^ with Ay = 0.187 in the rest frame V-band, 
and normmizing such that ^ P{Ay) = 1. 

3.6. Gravitational Lensing Degradation of the 
Luminosity Function 

As light from distant sources propagates to us, it is per¬ 
turbed by the numerous gravitational potentials of large 
scale structures along the way. The amplifying and de- 
amplifying the source as its photons are scattered into or 
out of our line of sight should result in a net zero effect, 
on average. However this lensing effect will add signifi¬ 
cantly to the dispersion in observed luminosity functions, 
and add an ever increasing disper sion with redsh i ft. W e 
use the simulated predictions from Holz & Linder (2005), 
which suggest the observed distribution in absolute mag¬ 
nitudes, given a gaussian-distributed internal dispersion 
in magnitude space, as a function of redshift will be 


am{z) = + (0.093 z)2, (7) 

where (Tint is the rest-frame dispersions given in Table 
We use this corrected dispersion in the adopted normally- 
distributed luminosity functions of Equation where 
PiMB) = N[MB,am{z)]. 


4. RESULTS 

The CCSNe rates from CANDELS and CLASH, and 
the weighted average of both, are shown in six redshift 
bins from z = O.lto 0 = 2.5in Table and Figure [5 
with Poisson errors based on the observed number o: 


events, derived from Gehrels (1986). The systematic er¬ 
ror range is shown in green, with two extinction scenarios 
as discussed below. For comparison, we show the rates 
from our D12 analysis, with nearly the same number of 
events per redshift bin as are in the CANDELS-fCLASH 
to z < 1.3. 


4.1. Systematic Errors in CANDELS+CLASH Rates 

Our approach to classification combines the tradition¬ 
ally systematic type and redshift uncertainties with the 
statistical error, making for a more complete, and per¬ 
haps fairer, assessment of overall rate uncertainty. How¬ 
ever there are other factors which clearly contribute to 
the systematic assumptions in our calculation. Our ad¬ 
ditional systematic error budget is shown in Table for 
each redshift bin, and includes the following factors: 


Type Eractions: While we have adopted the |Li et al^ 
(2011) volume limited CCSN type fractions tor com¬ 


pleten ess arguments, the fractions from Richardson et al. 
(2014) may also have a reasonably complete sampling of 
the CCSN diversity. Using the ASC-derived type frac¬ 
tions results in rates reduced by up to 20%. 


Peak Magnitudes: The value of the mean peak absolute 
magnitudes has a large impact on the discovery rate by 
shifting luminosity functions brighter or fainter against 
the survey efficiency function. Using the standard error 
in the mean peak absolut e magnitudes for each type from 


Richardson et al. (2014) leads to variances in the rates 
of as much as 45% in the highest redshift bins. 


K-correction: The k-corrections are determined from 
mean template SEDs provided by SNANA. We estimate 
the systematic uncertainty in these k-corrections using 
the dispersion in these mean corrections. It is important 
to note, the SNANA template libraries do not extend 
below ~3000A or beyond l/rm, making extrapolations 
from these templates to the synthetically observed pass- 
band somewhat uncertain. This is particularly worri¬ 
some in the lowest redshifts of this survey, at z < 0.4. 
It is difficult to estimate the amount of uncertainty in 
this spectral region, as there are surprisingly few NIR 
spectra with matched light curves for each main CCSN 
type in the literature. We use the error provided by 
the extrapolation of each template comprising the mean 
in the k-correction of each event (used to find U) over 
the passband widths. We find k-corrections add a 24% 
uncertainty to the (z) = 0.3 bin, and less than a 14% 
uncertainty in all higher-z bins. 
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TABLE 2 

Volumetric CCSN Rates From CANDELS and CLASH 



CANDELS 

CLASH 

CANDELS+CLASH 

Redshift 

Rate®' 

Ncc'" 

Rate 

Ncc 

Rate'^ 

Ncc 

0.3±0.2 


4.1 

o.oi_i 04 

5.7 


9.8 

0.7±0.2 

20 

14.8 

'^•^^-0.45 

0.9 

2-68ltol 

15.6 

1.1±0.2 

-1 -^o+l-OS 

YQ 

5.5 

-^•^y_0.78 

1.1 

-*-• ^^-0.71 

6.6 

1.5±0.2 

o qO + 1-98 

30 

6.1 

2 41+ 1^-3 

0.3 

q 9 c: + 2.03 

'^•^^-1.32 

6.4 

1.9±0.2 

q 1 C+3.24 

2.8 

q o-| +28.1 
'^•^-^-2.66 

0.2 

3.16t?'77 

3.0 

2.3±0.2 

c 1 c+6.59 
'^•^^-3.47 

2.6 

^■^-*^-5.30 

0.1 

Q 17+676 

-3.52 

2.7 


^ In units yr ^ Mpc ® 10 ^ h^Q with statistical errors. 

Erom Equation]^ in the intervals specified in the redshift column. 
Weighted averages of CCSNe rates. Uncertainties are standard errors in 
the weighted means. 


o 

a 





O 

Oh 




U 




Fig. 5. — Volumetric core-collapse supernova rates from CANDELS (blue squares) and CL ASH (yellow triangle s!. and the weight-average 
of both surveys (green circles). Points are offset horizontally for clarity. CCSNe rates from|Dahlen et al.|l|2012| white circles) are shown 
for comparison. Vertical error bars are the statistical uncertainties for each bin, while the horizontal error bars represent the width of the 
redshift bin. The green shaded region shows the systematic error range for the combined rates under two extinction assumptions: with 
variance in observed extinction, and under extreme limit scenarios discussed in Sectionj^ The right ordinate shows the scaled star-formation 
rate density, assuming i/i = Rcc/[^^^] with k = 0.0091 Mg 


Extinction in Host Galaxies: As stated in Section |3.5[ 
we consider the hosts of these supernovae to be largely 
normal galaxies, similar to the population of galaxies ob¬ 
served at low-z. The environmental extinction along the 
line of sight has similarly been assumed to be consis¬ 
tent with low-z CCSN observations, with Xy = 0.187. If 
we instead adopt an exponential distribution from more 
theoretical tr eatment of extin c tion, i.e., the MCMC pre¬ 
dictions from Riello & Patat (|2005 1 on the randomized 
extinction of events in spiral galaxies at different viewing 
angles (and setting Ay = 0.128), it would result in rates 
reduced by up to 20%, as shown in Table Our ap¬ 


proach to estimating the host extinction in the synthetic 
rest frame of the observed filter by scaling from the Ay 
using the Calzetti Law means that changes to the extinc¬ 
tion law, in Ry to within ±2, has a negligible effect on 
rate values. 

However, nature could be hiding a large fraction of 
CCSNe from us. By virtue of being tied to star-forming 
regions, CCSNe are also likely to occur in the dusti¬ 
est regions of galaxies. The strongest star-bursting 
galaxies, e.g. luminous and ultra-luminous IR galax¬ 
ies (U/LIRGS), which should produce the most CC¬ 
SNe, are also among the dustiest and UV/optically ob- 
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TABLE 3 

Systematic Error Budget in CANDELS+CLASH Rates^ 


Source 

(z> = 0.3 

{ z ) = 0.7 

( z > = 1.1 

(z> = 1.5 

{ z ) = 1.9 

(z) = 2.3 

Type Fractions 

+ 0.00 

- 0.34 

+ 0.00 

- 0.24 

+ 0.00 

- 0.35 

+ 0.00 

- 0.79 

+ 0.00 

- 0.76 

+ 0.00 

- 2.05 

Mb 

+ 0.64 

- 0.27 

+ 1.10 
- 0.26 

+ 0.81 

- 0.42 

+ 1.60 
- 1.02 

+ 1.63 
- 0.99 

+ 4.64 

- 2.65 

K-correction 

+ 0.71 

- 0.40 

+ 0.67 

- 0.00 

+ 0.16 

- 0.00 

+ 0.01 

- 0.01 

+ 0.11 

- 0.11 

+ 0.11 

- 0.11 

Extinction (Normal Galaxies) 

+ 0.00 

- 0.48 

+ 0.00 

- 0.36 

+ 0.00 

- 0.49 

+ 0.00 

- 1.06 

+ 0.00 

- 0.89 

+ 0.00 

- 2.02 

Extinction Limits 

+ 1.61 

- 0.54 

+ 3.14 

- 0.83 

+ 2.29 

- 0.60 

+ 4.42 

- 1.19 

+ 4.13 

- 1.35 

+ 8.31 

- 2.82 

Lensing Degradation 

+ 0.14 

- 0.00 

+ 0.33 

- 0.00 

+ 0.09 

- 0.00 

+ 0.03 

- 0.00 

+ 0.05 

- 0.00 

+ 0.29 

- 0.00 

ACDM 

+ 0.19 

- 0.18 

+ 0.16 

- 0.15 

+ 0.13 

- 0.12 

+ 0.24 

- 0.23 

+ 0.21 

- 0.20 

+ 0.49 

- 0.47 

Total Systematic (Normal Galaxies) 

+ 0.98 

- 0.78 

+ 1.26 
- 0.53 

+ 0.83 

- 0.74 

i + 1.62 

- 1.69 

+ 1.65 
- 1.55 

+ 4.67 

- 3.94 

Total Systematic (Extinction Limits) 

+ 1.89 

- 0.82 

+ 3.39 

- 0.92 

+ 2.43 

- 0.82 

+ 4.70 

- 1.77 

+ 4.45 

- 1.85 

+ 9.54 

- 4.40 


In units yr“^ Mpc“® /li^. 


scured ( Mannucd et al.|2007 ). This juxtaposition makes 
it difficult to account for expected dimming or loss of 
events, and is a point equally important to understand¬ 
ing the discrepancy between UV and IR-derived star for¬ 
mation rates and rate density histories. 

To consider a viable range, we evaluate two addi¬ 
tional extreme scenarios: a no extinction limit, where 
P{Av = 0) = 1 for all model contexts in tc', and a 
high extinction scenario, where a substantial fraction of 
high-z hosts are U/LIRG-like, obscuring large fractions 
of SNe from sight. For the hig h-extinction scenario, we 
use a prescription provided by Mattila et al. (2012) to 
account for the missing fractions by assessing SIN bud- 
get in U/LIRGs, and tracking the contribution of these 
Arp 299-like galaxies to the whole as a function of red- 
shift. We treat these high missing fractions as additional 
corrections to the rates derived for normal galaxies, us¬ 
ing the highest statistical bounds as the ‘high extinction 
limit’. These limits are meant as true extrema, not likely 
representative of environments of most GGSNe. How¬ 
ever, for completeness, we show these extrema in Figurej^ 
and Table [21 

Lensing Degradation: We have assumed an additional 
degradation in luminosity dispersions due to gravita¬ 
tional lensing. Assuming no lensing degradation, how¬ 
ever, results in a rates increased by less than 10%. 


instrumentation, detection thresholds, and other as¬ 
sumptions which go into their analysis. The reduced 
statistical uncertainty remains dominated by the sample 
uncertainty from this analysis, which could be further 
reduced with precise spectroscopic redshifts of the SN 
hosts, or a more consistent treatment of the D12 sample, 
in future studies. 

Also in Figure we compare our com bined 
GOODS-I-GANDE LS-I-GLASH rates to those from|Bot-| 


ticella et al.j (|2008l jz) = 0.21),[Gap pellaro et al. 
(z) = 0.20 and converted from SlNu), [Hazin et al. 


z) = 0.3), Melinder et al.| (2012 


iVlattila et a 


(2005 


(2DDg 

_ rrSo anc 

z) = 0.73), and Graur et al. )2011 (z) = 0.66), and 


to local z <0% rates from ^dappeharoet^ ([1^99), 


^igellaLetalJj 200g Smartt et al. 


15012), 'I'ayior et al. 120T3 ), anc 


Bot- 
,|Li e't al. 12011|) 


Graur 


et al. 1 2015). Assurnirig each is a valid measure of the 
GGbN rate at the mean or effective redshifts from each 
sample, whose certainty is characterized by statistical er¬ 
rors free from significant systematic offsets (i.e., ignoring 
systematic uncertainties), we determine a complete vol¬ 
umetric GGSN rate history, Rcc{z), using weighted av¬ 
erage rate measures in 5 equalized redshift bins in the 
range 0.21 < z < 2.34, with one additional weighted av¬ 
erage bin at z = 0.04 ± 0.04, as shown in Table [^ and 
Figure 


Cosmological Model: Throughout this analysis we have 
assumed the concordance cosmology, with = 0.27, 
Ha = 0.73, and Hn = 0.70. Using the range of uncer¬ 
tainty from Rest et al. (|2014 ), with 0.24 < Ha/ < 0.3 and 
0.70 < Ha < 0.76, results in as much as a 6% uncertainty 
in the resultant rates. 


4.2. Comparisons to Other CCSN Rates 

We combine the rates here with those from 
D12 to provide rates from the combined 
GOODS-I-GANDELS-I-GLASH programs (shown in 
Table and Figure [^, as they are similar in their 


4.3. Comparison to Cosmic Star Formation Rate 
Densities 

We compare our complete volumetric GGSN rate his¬ 
tory to cosmic star formation rate densities, SFR or 
ii^{z), in Figure]^ We show the co nsensus SFR derived 
from U V-f lR galaxy measures, from Madau fc Dickinson 


(2014) and Finkelstein et al. 
rived from IR luminous ga 


(|2014|), as well as a SFR 
axies (e.g., Ghary & Elbaz 
20011. The SFR, in units Mq yr“^ Mpc”"^, relates to the 
GGSN Rate, in units yr“^Mpc“^, by the fraction per 
unit mass of the initial mass function, IMF or ^{M), that 
is responsible for producing these GGSNe. We have de- 
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CJ 

a 






Fig. 6. — Rates from our group in comparison from other CCSN rates in the literature. Green circles: wei ghted average rates in six 
equalized redshift bins. Right Ordinate and Lines: star-formation rate density models, scaled to best match the jMadau &: Dickinson| SFR 
to all CCSN rate measures, with k = 0.0091 di 0.0017 . Also shown is the SFR model derived from the CCSNe rates (green) using the 
|Madau Sz Dickinson| parameterization. 


TABLE 4 

Comprehensive CCSN Rates 


Redshift Rate^ 


G O ODS+ GAND ELS-h GLASH: 


0.3±0.2 

9 1 q+0-80 

18.8 

0.7±0.2 

3.86t°;®® 

40.6 

1.1±0.2 

3.07lJ;0« 

17.6 

1.5±0.2 

Q 9c: + 2.03 

32 

6.4 

1.9±0.2 

3.16tf'77 

3.0 

2.3±0.2 


2.7 

RccizY-- 



0.04 ± 0.04 

n >79-1-0.06 
'^-0.06 


0.25 ±0.04 

1 oq-|-0.37 

i.OO_Q 29 


0.38 ±0.09 

1 81 ”l"0.31 


0.59 ±0.13 

Q Qi +0.95 


1.14 ± 0.42 

q 99 + 0.93 

'^•'^^-0.58 


1.93 ±0.37 

3.76t?;0i 



^ In units yr“^ Mpc“^ 10“^ h^Q 
with statistical eniprs. 

^ From Equationl2 in the intervals 
specified in the r^shift column. 

^ Weighted averages of CCSNe 
rates here and in the literature. 
Uncertainties are standard errors 
in the weighted means. 


termined this average scale empirically usi ng; a weighted 
least squares fit of the Madau & Dickinson (z) model 
to all CCSNe rate measures, Kcc{z), assuming they are 
related by 

Rcciz) = kh? -4)uviz). ( 8 ) 

With the scale, k, as the only free parameter, this com¬ 
parison results in fc = 0.0091 ± 0.0017 Mq^, where un¬ 
certainties are estimated from the reduced statical errors 
and the extent of the high and no extinction corrected 
rate values for the CANDELS-I-CLASH data only. We 
evaluate the goodness of fit from the reduced of the 
model fit to all r ate measures, and find xl = 1-4- A 
similar fit to the Chary & Elbaz (2001) 'ipiniz) model 
results in fc = 0.006 ± 0.002 Mq , but the goodness-of- 
fit is poorer, with xt = 5.3. A KS-test on the same 
model yields a D-statistic of 0.42, with a two-tailed p- 
value= 0.05, suggesting these distributions are unrelated 
to 95% confidence. 

In another comparison, we can attempt to predict the 
shape of tpiz) from the CCSN rates assuming the same 
parameterization as Madau & Dickinson (2014) for '0(z), 
such that 


ip{z) = 


A(l + z)^ 

{{l + z)/B)D + l 


( 9 ) 


We fit the above function to all rate data using a 
Levenberg-Marquardt least-squares algorithm, resulting 
in the green line and error region shown in Figure]^ with 
A = 0.015 ± 0.009, B = 1.5 ± 0.2, C = 5.0 ± 0.7, and 
D = 6.1 ±0.5. These values describe a SFR th at rises to, 
and declines from, z « 1 more steeply than the Madau fc 
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|Dickinson| model describes, and is likely to be inconsis- 
tent with SFR measures if extrapolated to very high-z. 

4.4. Comparison of Initial Mass Functions 
As shown in D12, from a numerical assessment of 


8 — 50 Mq progenitor stars with a Salpeter (19551 IMF, 
where: 

5OM0 

/ ^{M,z)dM 

_ ( 10 ) 

125M0 ’ 1 

/ M^{M,z)dM 

O.IMq 


k = 


one expects k = O.OO?)')^^^ Mq^, which is 25% lower 
than our fitted value but within the reasonable range of 
uncertainty, assuming the upper and lower mass bounds 
of CCSN progenitors remain in this range. 

There has been considerable concern on the absence of 
CCSN progenitors in nearby galaxies at masses > 20 Mq, 
and that, perhaps, this mass range represents a rough 
upper limit to the intrinsic nature of C CSNe progen- 

If so, 


itors (Smartt 
the expected 


2015 Eldridge et al. 2013|. If so, then 
raction of the liVll' responsible for CCSNe 


would drop to fc = 0.005 ± 0.002 (Salpeter). This 
value does not appear to be supported by the data (to 
about the 2cr level), but is difficult to completely reject 
given the various caveats in dust corrections for both the 
SFR and Rcc- There is also considerable uncertainty 
in the nature of the IMF, with the expectations (given 


8 — 50 Mq progenitors) of k = 0.0104 Mq^ for a Kroupa 
(2001) IMF, or fc = 0. 0079 fo r the steeper, high- 


mass (F = 1.45) IMF of Weisz et al. (2015), both equally 
consistent with the observed rates. 

As a spec ulativ e solut ion to the SFR discrepancy 
in galaxies, Dave (2008|) suggested a model in which the 
IMF evolves with redshift to increasingly favor the pro- 
ductio n of h igher mass stars in earlier cosmic epochs. 


Dave 


model, the 


In the 

which marks the change in the 


Kroupa 


characteristic mass, M, 
MF power-law slope for 


low-mass stars, moves to larger mass with redshift by 


M{z) = 0.5(1 


■ zf M, 


©• 


( 11 ) 


The model argues that tjj{z) models inferred from tracers 
of the light from predominantly high mass stars, assum¬ 
ing a present-day IMF, would overestimate the actual 
ip{z) at worsening rate with redshift. As the model sug¬ 
gests, scaling the observed ifuv{z) by an evolving frac¬ 
tion of high-mass contributors might reveal a truer 'if{z), 


one p ossibly traced by SN rates. We show the |Dave| 
(|2008|) V'(-z) niodel as the blue-dashed line in Figure]6| It 
is interesting to note, this implied 'if{z) would be similar 
in form (at least over the considered redshift range) to 
that expected fr om corrections by the evolu tion in stel¬ 
lar mass density ( Perez-Gonzalez et al.|2008 ), or perhaps 
from galaxy ‘downsizing’. 

Testing this tp^z) as we have for the previous mod¬ 
els, we find the best fit scaling to the Rcc{z) data 
with k = 0.0102Mq^. The Rcc{z) distribution is, 
however, inconsistent w ith th is ip{z) model, with a KS 
D = 0.58, rejecting the Dave| evolving IMF 'ip{z) model 
at the > 99.99% level. 

5. SUMMARY 

We have used a sample of approximately 44 CCSNe 
from the CANDELS and CLASH programs to extend 
volumetric CCSN rates to z = 2.5. Taken together with 
our previous results, we show a history of the CCSN rates 
from HST alone extending over the wide redshift range 
of 0.1 < z < 2.5, at a precision comparable to similar 
measures from the ground, at much lower z. In combi¬ 
nation with literature rates, we present a comprehensive 
CCSN rate history, with statistical uncertainties at or 
below systematic uncertainties. 

These composite rates are in good agreement with 
what would be expected from the consensus star- 
formation rate density history, assuming an average frac¬ 
tion of the IMF contributing to CCSNe progenitors of 
k = 0.0091±0.0017 Mq^ . The data are inconsistent with 
SFR histories derived solely from IR luminous galaxies, 
and less supportive of a progenitor mass range confined 
to the range 8 — 20 Mq, assuming a Salpeter IMF. The 
rate data also appear to rule out simple scenarios of an 
evolving (or relaxing) IMF. 

We present a model SFR histo ry derived from the 
CCSN rate data alone assuming a |Madau fc Dickinson| 
(2014|) p arameterization, whi ch is somewhat consistent 
with the Madau & Dickinson consensus ipiz) to z ~ 2.5, 
but is likely to have difficultly matching other ip{z) at 
even higher-z. 


This work is based on observations made with the 
NASA/ESA Hubble Space Telescope, delivered by the 
data archive team at the Space Telescope Science Insti¬ 
tute, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract 
NAS 5-265555. These observations are associated with 
programs GO-I2060 and GO-12099. 

HST 


APPENDIX 

CADENCE-AREAS OF THE CANDELS AND CLASH SURVEYS 

The full table of cadence-areas for each of the subfields in CANDELS, and all fields in CLASH, are shown in Tables 
and 1^ Table lists the areas of each template-search image pair within each target field that are co-added if they 
share the same cadence, or time between visits to the same pointin gs. Table [6] li s ts sim ilar combinations, but also 
combined for both parallel he lds of all 25 CLASH target helds. See |Grogin et al.| (|2011|) for a complete description 
of the CANDEL S fields, and Rodney et al. (2014) for a description ot how the CANDEL S SN survey varied from 
field to field. See Postman et aT 12012) tor a complete description of the CLASH fields, and Graur et al. (2014) for a 
description of how the CLASH SN survey varied from field to field. 
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TABLE 5 

Cadence-Areas eor the CANDELS MCT Fields 


Cadence 

(days) 

UDF 

GSA-Deep 

GSA-Wide 

GNA-Deep 

GNA-Wide SW 
(sq. arcmin) 

GNA-Wide NE 

EGS 

UDS 

GOS 

Total 

33 









4.60 

4.60 

45 


4.54 








4.54 

46 


44.08 

4.16 






13.78 

62.02 

47 


29.94 







78.12 

108.06 

48 


43.51 

4.42 


4.73 


4.24 


72.79 

129.69 

49 


89.95 

8.89 


4.49 


17.13 

89.39 

27.48 

237.33 

50 


65.13 

13.04 

15.05 

4.56 


63.13 

112.98 


273.89 

51 


34.31 


56.60 

3.27 


21.66 



115.84 

52 


106.88 


58.06 

11.39 





176.33 

53 


59.88 


36.11 

9.16 





105.15 

54 


24.87 


98.63 






123.50 

55 


44.07 


100.71 


38.66 




183.44 

56 


19.20 


78.08 

3.76 





101.04 

57 


5.00 


84.33 






89.33 

58 




24.66 






24.66 

59 




69.20 






69.20 

60 




56.43 






56.43 

61 


5.00 


10.15 






15.15 

62 








4.69 


4.69 

65 

5.11 



4.99 






10.10 

73 


4.40 








4.40 

85 



4.49 







4.49 

91 


24.94 








24.94 

92 


5.02 








5.02 

93 


34.84 








34.84 

95 


4.81 








4.81 

105 


4.89 








4.89 

136 


4.99 








4.99 

139 



4.43 







4.43 

140 




4.87 






4.87 

Sum= 

5.11 

660.25 

39.43 

697.87 

41.36 

38.66 

106.16 

207.06 

196.77 

1992.67 


TABLE 6 

Cadence-Areas for the 
CLASH MCT Parallel Fields 


Cadence ACS/WFC WFC3/IR 
(days) (sq. arcmin) 


6 

7 

8 
10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
30 

33 

34 

35 

39 

40 

41 
84 

Sum= 


4.81 

4.81 

4.81 

9.62 

76.99 

91.42 

48.12 

33.68 

14.44 

24.06 

9.62 

52.93 

43.31 

4.81 

4.81 

38.49 

28.87 

28.87 

19.25 

19.25 

19.25 

4.81 

4.81 

4.81 

19.25 

9.62 

4.81 

4.81 

634.92 


11.14 

11.14 

11.14 

22.28 

11.14 

189.35 

289.60 

89.11 

77.97 

44.55 

66.83 

44.55 

189.35 

155.94 

11.14 

100.24 

44.55 

66.83 

33.41 

55.69 

66.83 

11.14 


11.14 

22.28 


1648.72 
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